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1. Motivation and objectives 

Application of the method of LES to a turbulent flow consists of three separate 
steps. First, a filtering operation is performed on the Navier-Stokes equations to 
remove the small spatial scales. The resulting equations that describe the space- 
time evolution of the ‘large eddies’ contain the subgrid-scale (sgs) stress tensor that 
describes the effect of the unresolved small scales on the resolved scales. The second 
step is the replacement of the sgs stress tensor by some expression involving the 
large scales — this is the problem of ‘subgrid-scale modeling’. The final step is the 
numerical simulation of the resulting ‘closed’ equations for the large scale fields on 
a grid small enough to resolve the smallest of the large eddies, but still much larger 
than the fine scale structures at the Kolmogorov length. In dividing a turbulent 
flow field into ‘large’ and ‘small’ eddies, one presumes chat a cut-off length 6 can be 
sensibly chosen such that all fluctuations on a scale larger than 6 are ‘large eddies’ 
and the remainder constitute the ‘small scale’ fluctuations. Typically, 6 would be 
a length scale characterizing the smallest structures of interest in the flow. In an 
inhomogeneous flow, the ‘sensible choice’ for 6 may vary significantly over the flow 
domain. For example, in a wall bounded turbulent flow, most statistical averages 
of interest vary much more rapidly with position near the wall than far away from 
it. Further, there are dynamically important organized structures near the wall on 
a scale much smaller than the boundary layer thickness. Therefore, the minimum 
size of eddies that need to be resolved is smaller near the wall. In general, for the 
LES of inhomogeneous flows, we must consider the width of the filtering kernel 
6 to be a function of position. If a filtering operation with a nonuniform filter 
width is performed on the Navier-Stokes equations, one does not in general get 
the standard large eddy equations. The complication is caused by the fact that a 
filtering operation with a nonuniform filter width in general does not commute with 
the operation of differentiation. This is one of the issues that we have looked at in 
detail as it is basic to any attempt at applying LES to complex geometry flows. Our 
principal findings are summarized in this report. For details the reader is referred 
to Ghosal and Moin, 1993. 

In the field of subgrid-scale modeling, a large effort was invested in research 
on the dynamic localization model. The theoretical fundamentals have been set 
out in detail in our last report (Ghosal, Lund & Moin, 1993). In this report, we 
present some further results of tests performed on the model. These tests are on 
homogeneous turbulence, but tests on wall bounded flows have also been conducted 
by other members of the LES group and will be presented elsewhere (see the report 
by Cabot in this volume). 
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2. Accomplishments 

This section is divided into two parts. In §2.1, we summarize the theoretical work 
on the derivation of the basic LES equations for flows requiring nonuniform grids. 
In §2.2, some results from the ongoing tests of the dynamic localization model are 
presented. 

2.1 The basic equations for the LES of turbulent flows in complex geometry 

2.1.1 Nonuniform filtering in one space dimension (definition) 

Consider a field </>(£) defined in the domain (— oo,+oo). A filtering operation 
with a constant filter width A is defined by (Leonard, 1974) 

m = i /_ + °° G (^) (1) 

where G is any function with domain (— oc,+oo) and endowed with the following 
properties: 

(i) G(-0 = G(0 

(ii) /:~G(£)d£ = l 

(iii) G(£) — ► 0 as |£| — ► oo sufficiently fast so that all moments 

/ +oo 


( n > 0) exist. 

(iv) G(£) is very small (in some suitably defined sense) outside (— :|,+^). 
Some examples of possible filter functions are the ‘top-hat’ filter 


and the ‘Gaussian’ filter 


G(0 = 

U 

if Kl <> 

otherwise 

(2) 

G(0 = 

VI 

exp(-2£ 2 ). 

(3) 


(For a discussion of the various types of filters used in LES, see Aldama, 1990.) 

In situations where the domain might be finite or semi-infinite and a variable 
filter width is desirable, the definition (1) can have many possible generalizations. 
For example, a generalization of (1) when G is the top-hat filter might be 




i 

(A+(0 + A_(0) 


,«+A+(0 

/ <£(»?)<*»? 
J £— A_ (£) 


(4) 


where A+(£) and A_(f ) are positive functions and A+(£) + A_(f ) is the effective 
filter width at location ‘£\ For a finite or semi-infinite domain, A+(f ) and A_(f ) 
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Figure 1 . The space x, with variable filter width £(x), is shown mapped into a 
space £, with constant filter width A, by the ‘tan-hyperbolic map’ /(x) = tanh -1 x. 


must go to zero at the boundaries sufficiently rapidly so that (£ — A_(f ), £ + A+(f )) 
is always in the domain of <f>. 

A filtering operation with a nonuniform filter width (such as the one defined in 
(4)) does not in general commute with the operation of differentiation. It may be 
shown (Moin et al . , 1978), for example, that with the definition (4), 


d$_i( A + (0 + A_(0K i 

# ^ A+(0 + A_(() * (A + + A_) 


+ A+) 


d A+ 


Thus, 



+ ^-A.)^' 

(5) 

( 6 ) 


One would like to believe that the right-hand side of (5) would be small for some 
reasonable class of nonuniform filters, but this has never been conclusively demon- 
strated. This lack of commutativity between filtering and differentiation causes 
every spatial derivative operator in the Navier-Stokes equations to generate terms 
that cannot be expressed solely in terms of the filtered fields. Therefore, a ‘closure 
problem’ is introduced not only for the nonlinear terms, but for the linear terms 
as well. To remedy this situation we first propose an alternate definition for the 
filtering operation that is more general than (4). 

Let <f> be some field defined in a finite or infinite domain Any nonuniform 

grid in the domain a < x < b can be mapped to a uniform grid of spacing A in the 
domain [— 00 , + 00 ] by means of some mapping function 


f = /(*)■ 


Here f(x) is a monotonic differentiable function such that 


(7) 


f(a) = - 00 , 


( 8 ) 
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f(b) = +oo. (9) 

The nonuniform grid spacing 6(x) is clearly given by 

- m (10) 

(see figure 1). Clearly, if a (or b) is finite, (8)-(9) requires /' (a) (or /(b)) to be 
i nfini te so that 6(a) (or 6(b)) = 0. Since 6(x) can be regarded as a ‘local’ filter 
width (see equation (16)), the filtering kernel becomes a Dirac delta function at 
finite boundaries. This, of course, is an ideal limit that in a practical numerical 
computation can be achieved only approximately. Thus, our stipulation that the 
grid spacing be approximately equal to the filter width must break down when we 
are sufficiently close to a wall because in practice one cannot have an infinitely dense 
clustering of grid points at the boundary. 

The filtering operation is defined as follows. Given an arbitrary function */>(r), we 
first make a change of variables to f to obtain the new function </>(£) = <f-’( /~ 1 (<( ))• 
The function <?(£ ) is then filtered using the usual definition (1) appropriate for 
filtering on a uniform grid. Finally, we transform back to the variable r. Thus, 

$(x) = HO = £ jT G (^^) mdr), (11) 


or, on using (7), we have 

4>( x) = 1 ^*g( /( x) ~ /(v) ) 4iy)f'(y)dy. (12) 

Equation (11) or, equivalently, (12) is the definition we shall adopt for the filtering 
operation with a nonuniform filter width. For reasons that will become apparent in 
the next section, we will call this the Second Order Commuting Filter (SOCF). It 
should be noted that the definition (4) used by Moin et al is quite different from 
what one would get on substituting the expression (2) for the top-hat filter into 
(! 2 ). 

Example: In channel flow, one often uses the ‘tanhyperbolic grid’ (see for example 
Moin & Kim, 1982. The mapping function discussed below is slightly different 
from the version actually used in numerical computations since there one uses only 
a finite number of grid points whereas our filters become infinitely sharp at the 
walls). 

f(x) = tanh -1 x (13) 

where -hi > x > —1. (x = ±1 correspond to the channel walls.) From equation 
(12), the filtering operation is defined as 
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FIGURE 2a. The shape of the filter function G when G is a top-hat filter. 



FIGURE 2b. The shape of the filter function G when G is a Gaussian filter. 


where 

G(*,y) = ±G( m 2 na) )f'(v) (15) 

with f(x) = tanh 1 x. The function G(x,y) is plotted in figure 2A when G is a 
top-hat filter and in figure 2B when G is a Gaussian filter. If the approximations 
f( x ) “ f(y) » f'( x )( x ~ V) and f'(y) % f (x) for y near x are used in (15), we have 
on using (10) 

< i6 » 

Thus, as a first approximation, the filtering kernel G contracts in a self-similar 
manner on approaching a finite boundary (see figures 2 A and B). However, the exact 
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formula (15) has an additional higher order effect that causes the filter function G 
to become asymmetric near the wall giving more weight to points nearer the wall 
than further from it. The effect is most clearly seen in figure 2 A. 

2.1.2 Calculation of the commutation error 
Let us define the commutation error as 





(17) 


Then it is easily shown that 

CM = 1 /* c ( /(x) ^ — ) *'(*)/<») 


1 - 


/(*) 

f'(y) 


dy. 


(18) 


(The boundary terms are strictly zero due to the conditions (8) and (9).) 

It is convenient to introduce the new variable £ such that y is expressed implicitly 
in terms of ( through the equation 

f(y) =/(*) + A C (19) 


Equation (19) can be inverted by expressing y in a power series 


y = yo(C) + Ayi(C) + A 2 y 2 (C) + • • • 


( 20 ) 


where yo(C)> yi(C)> ' • • are functions to be determined. On substituting (20) in (19) 
and equating like powers of A, one obtains the expansion 




( 21 ) 


(Note: When the argument of any function is omitted, we imply that the function 
is evaluated at i x \ ) In terms of C, equation (18) may be written as 


/ + OG 

< 

-oc 


G(C)V> (y) 


l - 


/(*) 

f{y) 


d C 


( 22 ) 


where y is given by (21) (the limits of integration are obtained on using (8) and (9) 
in (19)). On expanding each of the factors in the integrand of (22) in Taylor series 
in A and collecting terms of the same order, we have, 


c[0] = C\ A + C 2 A 2 + • • • — C|6 + C2^ 2 + 


where 


ci = / <G(C)dC, 

f J -OG 


(23) 

(24) 
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c 2 


2/7V +//*>*- 3/" V / +oc 


2 /' 


>4 


/: 


C 2 G(C)dC, 


(25) 


Cj — ci/ and c 2 = c 2 (/ ) 2 . Since G(() is symmetric, ci = c, = 0. Thus, the 
commutation error C[t/>] ~ 0(6 2 ). 

In an LES, the grid spacing is approximately equal to the ‘filter-width’, 6. If a 
second order numerical scheme is used to represent the derivatives, the finite differ- 
encing error is then of the same order as the error due to the lack of commutativity 
of the differentiation and the filtering operations. Therefore, in an LES of an in- 
homogeneous turbulent flow using a second order finite differencing scheme, the 
filtering operation can be considered to commute with the differentiation operation 
to within the accuracy of the numerical approximation. This is our most important 
result from the point of view of practical application. 


2.1.3 Spectral distribution of the commutation error 
Let us substitute 

= 4'k exp(ikx) (26) 

and (21) in (22). Then, 


C[*!>] = ik%l> f G{ C) 
J — oo 


1 - 


/(*) 


r(* + f + -)J 

On expanding the integrand of (27) in a power series in A, we have 


exp dC (27) 


cub) 

“p = Fo(kA) + A^i(l'A) -| (28) 

where To, T\, . . . contain A only in the combination kA. 

Now A < 1, but kA may be as large as order one. Thus, we may neglect the 
successive terms in (28) in comparison to the first term to obtain 


gM 

*/> 


~ To(kA). 


On evaluating To(kA) from (27), we have 


m 



(29) 


(30) 


where T is defined by 



+ oo 


CG(C)exp(ixC)dC. 


(31) 
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Equation (30) can also be written as 

m 


xj) 


= -ikS 


7{k6) 


(32) 


where 6 is the local filter width as defined in (10). On expanding the exponential in 
(31) in a Taylor series, we see by virtue of G{Q ) being symmetric that J-(kS) ~ kS 
so that |C[V’]/V’I ~ (l'<5) 2 as shown in the last section. 

Comparison with finite differencing errors is facilitated if the commutation error 
is expressed as a ‘modified wave-number’. If C\tp] were zero, we will have for the 
function (26) 


ax ax 


Therefore, if we define a ‘modified wave-number’ k by 


^ = t 
dx 


(33) 


(34) 


then the departure of k' from k is a measure of the commutation error. On making 
the change of variable (19) in (12), we obtain 


t Kx)= ( + °° G(04>(y)dC. 

J — OO 

On substituting (26) in (35) and on using the definition (34), we obtain 


(35) 


k J-k G(( ) exp( iky(x, ( ))d( 


(36) 


where y has been expressed as a function of £ by inverting (19) for each fixed value 
of x. Equation (36) is an exact result. A simplified asymptotic form is obtained 
on substituting the expansion (21) in (36) and dropping all terms such as A(fcA), 
A 2 (fcA) • • • which are negligible compared to kA. Thus, 


f , A r /r<G(oex P ^ 

= 1 - A 


dC, 


k ' f' 1 /rG(C)exp(^ 

Since G(£ ) is a symmetric function, (37) simplifies to 

k 


d C 


(37) 


f" 


K 

f' 2 C)cos( 

:?) 

\dC 


( 38 ) 
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0 2 4 6 6 

wavenumber, k 


FIGURE 3a. The modified wavenumber (k ) for the commutation error (— ) 
compared to that of the centred differencing error (— — — ). The symbols k X' are 
the result of using the approximate formula (42). Here A = 27r/16, x = 0 (channel 
center), and the maximum wavenumber is 7 r/6 where 6 is the local filter width. 


Example : Let us consider the top-hat filter defined by (2) together with the tan- 
hyperbolic map f(x) = tanh” 1 x. For this map, equation (19) can be inverted to 


give 


x + tanh 
1 + x tanh A £ 


(39) 


On substituting (39) in (36), we obtain 



P(ih)G(0 


I — tanh 2 AC 
( 1 + jt tanh AO 2 


/^exp(/A-j/)G(CVC 


( 40 ) 
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wavenumber, k 

FIGURE 3b. Same as figure 3A for x = 0.95 (close to channel wall). 

with y given by (39). As a typical example, we consider a channel whose walls are at 
x = 0 and x — 2n with 16 grid points in the spanwise direction. Thus, A = 2 tt/ 16. 
The integral (40) can be evaluated numerically. The result is shown in figure 3. The 
modified wavenumber for the second order central difference scheme is given by 


k sin (k6) 
~k = k6 


(41) 


Equation (41) is also plotted in figure 3 for comparison. The asymptotic formula 
(38) can be evaluated analytically in the case of the tan-hyperbolic map and top-hat 
filter. A straightforward computation gives 
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Equation (42) is also shown in figure 3. The agreement of the asymptotic result 
(42) with the exact result (40) is seen to be very good. 

The asymptotic formula (38) can be written in the following convenient form: 

k R = k, (43) 


*/ = \jj F(kS) 

where 6 is the local filter width given by (10), 


F(x) — x 


/_ ”CG(C) sin (gC)<*C 

G(Qcos(xC)dC 


and the suffixes R and I denote real and imaginary parts respectively. Thus, to a 
very good approximation, the commutation error is seen to be purely dissipative 
in nature in contrast to the central differencing error, which is dispersive. The 
commutation error vanishes in regions where 6 — 0 (such as at the center of a 
channel). From the example above, it is clear that for the top-hat filter, 


, k6 ( k6\ (k6) 2 ( kS ) 4 

fW-l- T «^ T ) = LJ- + L-i- + .... (46 ) 

A simple calculation shows that for the Gaussian filter (3), 

W = (47) 

2.1.4 Higher order corrections to the commutation error 

We have shown in the previous section that C[4'\ ~ 4 y Fo(kA) and Tq ( kA) ~ 
(A:A) 2 at leading order in kA. In this section, we shall attempt to approximate 
the commutation error C[ift] by an expression involving 4? and its derivatives such 
that the residu al is of order (&A) 4 . The procedure can be readily generalized to 
represent ij>'(x) in terms of (x) and higher derivatives of 4’(x) such that the error 
in the approximation is at most of order (&A) 2m w r here m is any positive integer. 

We have, on expanding the exponential in (31) and noting that G(Q is a sym- 
metric function, 


^ C 2 G(CK - l ±.j c* 


G(CK + 


Substitution of (48) in (30) gives 


(kA)1 Uy*L 


eG(C)dC + 0(kA)\ 


(49) 
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From (35), 


V>(i) 




(50) 


On substituting ip = ip k exp(ikx) in (50) and differentiating twice with respect to 
x, we have 

r / * A T 

(51) 


f '(x) = -Pip [l - + 0(*A) 4 

2/ 


where 


/ +oo 

I 

-oo 


C 2 G(CMC 


and A < 1 has been assumed. Equation (51) implies that 

-k 2 ip = f(x) + 0(k&) 2 . 


(52) 


(53) 


On substituting (53) in (49) we get 


C[iP} = aj^A 2 f(x) + 0(kA) i . 


(54) 


Thus, 


dip _ dtp f 

dx dx a f' 3 dx 2 


+ O(itA) 4 . 


(55) 


The procedure can be continued to extend the accuracy of the representation to any 
order in fcA. Equation (55) can also be written in terms of the local grid spacing 
S(x) as follows: 


difr 

dx 


- it - 

dx 


A 


(5 


) 


dx 2 


+ 0(k6)*. 


(56) 


Equation (56) was established only for the function (26). However, it is clearly valid 
for any linear superposition of functions of the type (26), that is, any function that 
admits a Fourier representation. 


2.1.4 Generalization to three space dimensions 

These results can be generalized to three space dimensions (see Ghosal and Moin, 
1993 for details). Let us consider curvilinear grids defined by the co-ordinate sur- 
faces 

H(x) = constant (57) 


where x are rectilinear co-ordinates in physical space. Let us also introduce a new 
1 computational space’ X through the map 


X = H(x) 


(58) 
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0.5 1.0 1.5 2j0 2.5 3.0 

time 


FIGURE 4A. Decay of energy in a simulation of the Comte- Bellot and Corrsin 

experiment with 48 3 grid points. Large scale energy ( ); sgs energy (k) ( ), 

as predicted by DLM(G). This is compared with an identical computation with the 

DLM(C) (+), the DM ( ), and with the model turned off ( ). Experimental: 

resolved energies (•); subgrid energies (■) 

which maps the physical space domain into R 3 meshed with a uniform grid of 
spacing A. The filtering operation is now defined as 

*x> = ^ *(*:> ) tfww (59) 

where J(x) is the Jacobian of the transformation (58). 

It can be shown that for the general three dimensional filter (59), equation (54) 
takes the form 

C k [ti>i = — oA 2 r + o(|k|A) 4 (60) 

where 

r kmp = fc mJf (H(x))A, tf (H(x))jy il *(x) (61) 

and h is the inverse of H. Further, an expression for Ck[4>] in terms of 0 can be 
written down to any order in A. Thus, 

dki J = (dk — otA 2 rkmnC% nn H W 1 - ( 62 ) 

2.2 Tests of the dynamic localization model 

Two versions of the dynamic localization model (DLM) were introduced in last 
year’s Annual Research Briefs (Ghosal, Lund & Moin, 1993). The first of these 
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wavenumber 

FIGURE 4B. The spectra at three times for the same computation as in figure 4A. 

DLM(G) ( ); DM (— - -); DLM(C) (+); no model ( ). The solid circles, 

squares, and triangles are the experimental values at the three experimental points 
in figure 4A. The initial conditions are chosen so that the spectra at the earliest 
time matches the experiment. 



time 

FIGURE 5A. Time dependence of energy in a forced simulation of homogeneous 
isotropic turbulence at infinite Reynolds number on a 32 3 grid using the DLM(G). 

Energy per unit volume in the resolved scales ( ); subgrid- scale energy (— ); 

Total energy ( ). 

(which we will call the ‘constrained’ dynamic localization model or DLM(C)) con- 
strains the Smagorinsky coefficient to be nonnegative and therefore does not exhibit 
the phenomenon of backscatter - the (locally) reverse transfer of energy from the 
subgrid to supergrid scales. The second formulation (which we will call the ‘gen- 
eral’ dynamic localization model or DLM(G)) removes this restriction. However, 
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FIGURE 5B. Equilibrium spectrum in a forced simulation of homogeneous isotropic 
turbulence at infinite Reynolds number on a 32 3 grid. Kolmogorov’s five-thirds 
law would correspond to a horizontal straight line with an ordinate equal to the 
Kolmogorov constant C K % 1.5. DLM(G) ( -); DLM(C)(- - -); DM ( ). 
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a, 



time 


FIGURE 5C . Fraction of backseat ter in a forced simulation of homogeneous isotropic 
turbulence at infinite Reynolds number on a 32 3 grid using DLM(G). Volume av- 
eraged rate of energy transfer from subgrid to supergrid scales as a fraction of the 

net transfer ( ); Fraction of points experiencing backscatter at a given instant 

of time (— — — ). 


the generality is gained at the cost of added computational expense. The so called 
‘volume-averaged’ version of the dynamic model which can be rigorously derived 
as a special case of the variational formulation (Ghosal, Lund & Moin, 1993) will 
be simply referred to as the dynamic model (DM). It is useful only for flows with 
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homogeneous directions. 

Previous tests had shown (Ghosal, Lund k Moin, 1993) that while DLM(C) and 
DM used in a 32 3 LES simulation of the Comte-Bellot and Corrsin experiment 
(Comte-Bellot and Corrsin, 1971) gave excellent predictions, the results obtained 
with DLM(G) were not as satisfactory. This discrepancy was more thoroughly 
investigated. It was found that at a resolution of 32 3 , the subgrid-scale energy (k) 
was almost equal to the total energy in the LES field. As such a situation is not 
fully consistent with the basic assumption of LES that the resolved scales carry 
most of the energy, the simulation was repeated on a more refined (48 3 ) mesh. 

At this resolution the total resolved energy is about thrice the subgrid-scale en- 
ergy, and the performance of the DLM(G) is greatly improved. The performance 
of DLM(C), DLM(G,) and DM are compared in figure 4 starting with identical ini- 
tial conditions consistent with the experimental data at the initial time. All three 
models are seen to give comparable results in good agreement with the experiment. 
The result of running the simulation with the model turned off is also shown in the 
same figure. The crucial role played by the model in these simulations is apparent. 

This was followed up with a simulation of forced homogeneous isotropic turbu- 
lence at infinite Reynolds number using a 32 3 grid. The forcing was provided by 
the simple artifice of resetting the six modes corresponding to k = (±1,±1,±1) 
(k is the wavenumber vector) to fixed values at the end of each time-step. The 
achievement of steady state was monitored by plotting the energy as a function of 
time and also by looking for time independence in the shape of the spectra. Fig- 
ure 5 A illustrates the extent to which a ‘steady-state’ was achieved during a run 
using the DLM(G). Though a plateau was reached in the total energy, some oscil- 
lations remained. These could probably be reduced by running the simulation for 
a longer time, but part of it might be an artefact of the forcing scheme. If a strict 
Kolmogorov spectrum was achieved, the sgs energy would be 

roc o 

k — I C K e 2/i K~ 5/3 dK = -C s f 2/5 «; 2/3 (63) 

where t is the dissipation rate, K m is the wavenumber magnitude at cutoff, and 
C K 1.5 is the Kolmogorov constant. In our simulation, n m = 14 (the number of 
full energy shells) and e « 0.07. If these numbers are substituted in (63,) we find 
k « 0.03. It is reassuring to note that the computed sgs energy (see figure 5A) is 
within a factor of two of this figure. (Of course, an accurate quantitative prediction 
of subgrid-scale energy cannot be expected from a large eddy simulation.) 

The quantity e -2 / 3 /c 5 / 3 E(k), where E(k) is the energy spectrum and e is the sgs 
dissipation, is plotted in figure 5B for runs using each of the three models. No time 
averaging has been performed, so there is some sampling fluctuations especially at 
the low wavenumbers. All three models show good agreement with Kolmogorov s 
law with a Kolmogorov constant close to 1.5. The rising ‘tail’ in the DLM(G) 
simulation at high wavenumber appears to be real since the high wavenumber part 
shows negligible fluctuations among different samples. On the other hand, the 
DLM(G) appears to give a better plateau at intermediate wavenumbers at a value 
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very close to C K cs 1.5. Again, this behavior was seen when the simulation was 
repeated with a different initial spectrum, so the effect is probably real and not just 
an accidental sampling fluctuation. 

The percentage of ‘backscatter’ is shown in figure 5C using two different measures. 
The solid line is the rate of transfer of energy from the subgrid to supergrid scales 
(integrated over the volume of the box) as a fraction of the total transfer. The 
dashed line is the fraction of points engaged in this ‘backscatter’ of energy at a 
given time. The results on backscatter can be tested using a DNS database of 
a high Reynolds number simulation of homogeneous turbulence, but this has not 
been done yet. However, Piomelli et al (1991) found, using a DNS of channel 
flow at Re = 3300, that the fraction of backscattered energy is about 10-20 percent 
while the fraction of backseat tering points is about 20-30 percent. These numbers 
depend on the type of filter used. The figures quoted here are for the ‘box filter’ 
(see figure 9 in the paper of Piomelli et al.) which is what we are using as our 
‘test-filter’. Figure 5C is roughly consistent with Piomelli et aV s results. These 
simulations were repeated using a more refined (48 I * 3 ) mesh. None of the above 
results were changed in any significant way. 

If one has to choose a model for an LES of a complex flow, then, based solely on 
these tests, the DLM(C) will perhaps be the model of choice since it is the simpler 
of the two available local models and gives results that are just as good as the more 
elaborate DLM(G). The extra effort invested in modeling backscatter does not seem 
to significantly improve predictions of quantities of practical interest. However, in 
more complex flow situations where backscatter is known to play an important role, 
the DLM(G) might turn out to be the preferred model. This issue will be settled 
in the near future as the results of applying these theories to complex flows become 
available. 

3. Future plans 

It is expected that these ideas will be useful in the near future in the context of 
the application of the dynamic localization model to LES of complex flows. 

An explicit knowledge of the nature of the filtering operation was not needed in 
the past in sgs modeling. However, all models in the ‘dynamic model’ category rely 
on an explicit ‘test’ filtering operation (which is supposed to be self-similar to the 
grid-level filter) to extract useful information about the subgrid scales. Thus, the 
issue of the commutativity of the filtering and differentiation operations is important 
not only for providing a firm theoretical basis for LES in complex geometry, but 
also for practical calculations which make use of the dynamic model. 

The DLM should be applied to progressively more challenging flows to evaluate its 
performance. Faster and more efficient algorithms for solving the relevant integral 
equations need to be examined. 
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